Constraining relativistic models through heavy ion collisions 
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Relativistic models can be successfully applied to the description of compact star properties in 
nuclear astrophysics as well as to nuclear matter and finite nuclei properties, these studies taking 
place at low and moderate temperatures. Nevertheless, all results are model dependent and so far 
it is unclear whether some of them should be discarded. Moreover, in the regime of hot hadronic 
matter very few calculations exist using these relativistic models, in particular when applied to 
particle yields in heavy ion collisions. In the present work we comment on the known constraints 
that can help the selection of adequate models in this regime and investigate the main differences that 
arise when the particle production during a Au+Au collision at RHIC is calculated with different 
models. 
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I. INTRODUCTION 

Relativistic models have been widely used in order to 
describe nuclear matter, finite nuclei and stellar matter 
properties. Many variations of the well known quantum 
hadrodynamic model [l| have been developed and used 
along the last decades. Some of them rely on density 
dependent couplings between the baryons and the mesons 
[2, H, 0, IS @1 while others use constant couplings @, H, @ • 
Still another possibility of including density dependence 
on the lagrangian density is through derivative couplings 
among mesons and baryons [lCj, LLlI, [l2| or the coupling 
of the mediator mesons among themselves [l3[ . 

All these quantum hadronic models have under control 
the fitting of the infinite nuclear matter binding energy at 
the experimental saturation density (po). Their results 
are around 16MeV and 0.15/m -3 , respectively. How- 
ever at the level of nuclear matter, they predict different 
values for other physical quantites, as for instante, incom- 
pressibility (K), effective nucleon mass, scalar and vector 
potentials, critical temperature (T c ), etc. If applied for 
finite nuclei, these models also present discrepancies on 
nuclear spectra. The efforts to find some correlations 
among observables of infinite nuclear matter and finite 
nuclear results, lead to interesting conjectures. One is 
that there is a correlation between the effective nucleon 
mass (M*) at the saturation nuclear matter density and 
the spin-orbit splitting for several nuclei 14j , for different 
quantum hadronic models. This correlation states that 
good theoretical L-S splittings for several finite nuclei are 
obtained if M* lies between 0.58 and 0.62. Other corre- 
lation was also proposed to connect T c with K, M* and 
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There is a strong correlation between compact star 
properties and some of the nuclei properties. Relativis- 



tic model couplings are adjusted in order to fit expected 
nuclei properties such as binding energy, saturation den- 
sity, compressibility and energy symmetry at saturation 
density, particle energy levels, etc. Once the same rel- 
ativistic model is extrapolated to higher densities as in 
stellar matter or higher temperatures as in heavy-ion col- 
lisions or even to lower densities as in the nuclear matter 
liquid gas phase transitions (lq . they can and indeed 
provide different information. Hence, experimental con- 
traints obtained either from polarized electron scattered 
from a heavy target, from heavy-ion collisions at differ- 
ent energies or from astronomical observations are very 
important in order that adequate models are chosen and 
inadequated ones are ruled out. In what follows we dis- 
cuss some of the already existing constraints and the re- 
lation between equations of state (EoS) used to describe 
stellar matter and hadronic matter. 

The relation between neutron star properties which are 
obtained from specific EoS and the neutron skin thickness 
has long been a topic of investigation in the literature. 
The difference between the neutron and the proton radii, 
the neutron skin thickness, is linearly correlated with the 
pressure of neutron matter at sub-nuclear densities. This 
is so because the properties of neutron stars are obtained 
from appropriate EoS (very isospin asymmetric due to 
the (3- equilibrium constraint) whose symmetry energy 
depends on the density and also controls the size of the 
neutron skin thickness in heavy and asymmetric nuclei, 
as 208 Pb, for instance. In [ItJ it was shown that the mod- 
els that yield smaller neutron skins in heavy nuclei tend 
to y ield smaller neutron star radii due to a softer EoS. In 
18] it was shown that the neutron skin thickness indeed 
give hints on the equations of state that are suitable to 
describe neutron stars. A correlation between the slope 
of the symmetry energy and the neutron skin thickness, 
previously found for Skyrme-type models [19|, was also 
observed within relativistic models [18j • Unfortunately a 
precise measurement of the neutron skin thickness is still 
under way. 
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Neutron stars are believed to have a solid crust formed 
by nonuniform neutron rich matter in /3-equilibrium 
above a liquid mantle. In the inner crust nuclei coex- 
ist with a gas of neutrons which have dripped out. The 
properties of this crust as, for instance, its thickness and 
pressure at the crust-core interface depend a lot on the 
density dependence of the EoS used to describe it [20| . 
On the other hand, it is well known [2l|, |22| that the ex- 
istence of phase transitions from liquid to gas phases in 
asymmetric nuclear matter (ANM) is intrinsically related 
with the instability regions which are limited by the spin- 
odals. Liquid-gas phase transition in ANM can lead to an 
isospin distillation phenomenon, characterized by a larger 
proton fraction in the liquid phase than in the gas phase. 
In [23| the spinodal sections for two relativistic models, 
one with constant couplings known as NL3 [7J and an- 
other with density dependent couplings we usually refer 
to as TW, were obtained. The curve that represent the 
proton versus neutron densities for /3-equilibrium matter 
was represented in the same plot and one can see that it 
always crosses the spinodal section at T = 0. This means 
that a liquid-gas phase transition occurs at the crust of 
a cold neutron star giving rise to a nonhomogeneous re- 
gion. For higher temperatures this effect generally does 
not occur meaning that the outer layer of a hot compact 
star is homogeneous. For T = 10 MeV, for instance, the 
two models investigated gave contradictory information 
as for the NL3 the crossing was observed and for the TW 
it was not. More information is then required so that the 
correct picture can be depicted and the wrong EoS ruled 
out. 

Compact star properties can be computed from the so- 
lution of the Tolman-Oppenheimer-Volkoff (TOV) equa- 
tions [2~I1 ] once they are supposed to be spherically sym- 
metric and static. The input to the TOV equations 
is the EoS chosen to describe the stellar matter. The 
determination of neutron star properties obtained from 
the measureament of the gravitational redshift of spec- 
tral lines produced in neutron star photosphere provides 
a direct constraint on the mass-to-radius ratio. Some 
time ago a redshift of 0.35 from three different tran- 
sitions of the spectra of the X-ray binary EXO0748- 
676 was obtained in (25|. This redshift corresponds to 
M/R = 0.15M e /Km. The IE 1207.4-5209 neutron star, 
which is in the center of the supernova remnant PKS 
1209-51/52 was also observed and two absorption fea- 
tures in the source spectrum were detected [26j. These 
features were associated with atomic transitions of once- 
ionized helium in the neutron star atmosphere with a 
strong magnetic field. This interpretation leads to a read- 
shift of the order of 0.12-0.23, considerably lower than the 
one in 25]. This rcadshift imposes another constraint to 
the mass to radius ratio given by M/R — 0.069Mq/ Km 
to M/R = O.II5M0/ Km. However, the interpretation of 
the absorption features as atomic transition lines i n 1 261] 
is controversial: an alternative interpretation [27|, |28|] is 
that the absorption features are cyclotron lines, which 
imply no obvious constraint on the EoS. In previous 



works [29|, |30j, l3l| many EoS were tested against these 
two possible constraints. While all EoS studied so far 
are consistent with the measurement proposed by I2q| . 
many of them are excluded if the redshift found in [25j 
provides the only possible constraint. As an example we 
plot the mass and radii obtained from different EoS in 
Fig. [TJ where we have added the line corresponding to 
the more restrictive constraint (top straight line). 
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FIG. 1: Mass-radius plots for one hadronic, three hybrid and 
two quark stars. 

The hadronic star [32j w &s built with the inclusion 
of 8 baryons within the non-linear Walecka model and 
parametrization GM1 0] given in fl] The hybrid stars 
were obtained with the same model and parametrization 
as above for the hadron phase and different models for 
the quark phase. Hybrid 1 [33j mixed phase contained 
the MIT bag model f|[ with Bag = (160 MeVVV hybrid 
2 33] the color flavored locked phase (CFL) [36| again 
with Bag = (160 MeV) 4 and hybrid 3 H3 the Namb- 
Jona-Lasinio (NJL) model [Hj]. The quark star proper- 
ties were obtained for the CFL with the same bag pa- 
rameter as above and called quark 1 and with the NJL 
model and called quark 2 [33[ ■ One can notice that while 
hybrid 2 star is excluded, hybrid 1 and hybrid 3 stars 
remain valid. The hadronic and quark stars remain as 
possible candidates. 

Electron-positron pairs can be emitted from bare quark 
stars 34] and this dominant emission requires the exis- 
tence of a surface layer of electrons tied to the star by a 
strong electric field. While the electron chemical poten- 
tial of a quark star described by the MIT bag model is 
very low (less than 20 MeV), the NJL model gives much 
higher values reaching 100 MeV inside the star [39J]. 

A very recent paper reported a series of pos- 
sible constraints used separately to analyse relativistic 
and non-relativistic equations of state applied to neu- 
tron stars. The authors mention a relation derived in 
[4ll ] that correlates gravitational and baryonic neutron 
star masses and the measurement of Pulsar A moment 
of inertia [12] which excludes some of the standard mean 
field relativistic models. 

From the above considerations, one can see that differ- 
ent relativistic models in the description of compact stars 
provide different information. In order to identify the 
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possible constituents of compact objects and, in this way, 
rule out some of the relativistic models, more astronom- 
ical observations are necessary. 

Another attempt in the direction of choosing appropri- 
ate models is discussed in the present work in the high 
temperature regime produced in heavy ion collisions. For 
this reason we consider Au-Au collisions at RHIC/BNL 
and analyse the hadron abundances and particle ratios in 
order to determine the temperature and baryonic chem- 
ical potential of the possibly present phase transition. 
The possibility that a quark-gluon plasma (QGP) could 
be formed in heavy ion collisions arose when quantum 
chromodynamics (QCD) at finite temperature and high 
densities became a topic of increasing interest due to the 
discovery of asymptotic freedom about 30 years ago. In 
cosmology, the relevant conditions for QGP formation 
occur 10 fxs after the Big-Bang theory; nuclear matter 
first appears after about 1 ms. In laboratory searches for 
QGP, in large colliders around the world (RHIC/BNL, 
ALICE/CERN, GSI, etc), experimentalists are trying to 
do the opposite: to convert hadronic matter at suffi- 
ciently high temperatures into QGP. 

A first step towards the interpretation of data collected 
during a Pb+Pb collision at the SPS was taken in [43| . 
where the authors have used a statistical model which 
assumes chemical equilibration to find the temperature 
and baryon chemical potential that provide a best fit to 
the data obtained by the NA49 [HJ and WA97 [H| co- 
laborations. They concluded that the chemical freeze- 
out temperature was T = 168 ± 2.4 MeV and the cor- 
responding baryon chemical potential was ps = 266 ± 5 
MeV. Lattice Monte Carlo simulations of QCD at van- 
ishing baryon density gives a very similar estimate, i.e., 
T = 170 ± 8 MeV Latter on, the same model was 
used to analyse the particle production yields measured 
in central Au+Au collision at RHIC [47j and the au- 
thors obtained T = 174 ± 7 MeV and p B = 46 ± 5 MeV. 
In both cases 0, |47} they used the free Fermi and Bo- 
son gas approximations where the interaction among the 
baryons and mesons were neglected and there were only 
the above mentioned independent parameters, i.e., the 
temperature and baryon chemical potential. Three con- 
servation laws were imposed: of the baryon number, the 
strangeness and the electric charge. While the temper- 
atures found were of the same order, the very different 
chemical potentials imply different baryon density. In a 
more recent paper [48| , a simple relativistic selfconsistent 
chiral model, where the nuclear interaction was included, 
was employed to analyse the Au+Au collision at RHIC 
once again. The results depend on the parametrizations 
used but they found a lower temperature around T = 155 
MeV and the baryon chemical potential of the order of 
Hb = 51 MeV. The authors claim that the fitted chemical 
freeze-out temperatures and chemical potentials depend 
on the order of the phase transition and suggest that at 
RHIC the system emerges after the phase transition. 

In the present work we revisit the same data and anal- 
yse them with the help of more sofisticated relativis- 



tic models within a mean field approach. We use four 
parametrizations of the non-linear Walecka model [l|, 
namely NL3 0, TM1 [|, GM1 and GM3 0], one model 
with implicit density dependence through meson field 
couplings, the NLwp [l3T ] and two different parametriza- 
tions of a density dependent hadronic model, the TW 
[3j and the DDME1 [g]. Similar investigations were per- 
formed in [49[ | within one of the parametrizations men- 
tioned below (GM1) and the best fit resulted in T = 143.9 
MeV and ps = 25.07 MeV. We need to understand how 
model dependent these numbers are. Hence, for com- 
parison, we also investigate the results obtained within 
a model containing the eight lightest baryons (baryonic 
octet) without any interaction and then with the bary- 
onic octet plus extra ten baryons (baryonic decuplet) also 
without interaction to try to understand how important 
the interaction at low densities are. We implemented a fit 
based on the minimum value of the quadratic deviation 
as in [43] in order to obtain the temperature and chem- 
ical potential for each model. The paper is organized as 
follows: in section II we show the lagrangian densities of 
the models considered and describe the formalism used; 
in section III we present and discuss the results; in section 
IV we draw our final conclusions. 



II. FORMALISM 

In the following we review the basic formulae of the 
models we use. Just the Lagrangian density and the 
main expressions for densities and chemical potentials 
are given. If one wants to follow the analytical calcula- 
tions in detail, please refer to [E EM US HH j among many 
other papers in the literature. 



A. Different parametrizations for the non-linear 
Walecka model 

The Lagrangian density that incorporates many 
parametrizations of the non-linear Walecka model 
(NLWM) [J, 0, H, II and also the extra non-linear a - p 
and ui — p couplings fl3j reads 



£ = ^ B b» C* 5 " - 9vbV- - 9pB t ■ 

B 



-eii±M#|-(M- M 
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W[k s g 2 sB 4> 2 + A v g 2 vB V p V% 



(1) 



where B represents the lightest 18 baryons, cf>, V^, 
b M and N 1 are the scalar-isoscalar, vector-isoscalar and 
vector-isovector meson fields and the photon field respec- 
tively, ty^, = d^Vv-dyVp , B M „ = d^v-dvh^-TpiJopy, 
by), = d^Ay — <9„A M and T3B is the third component 
of the baryon isospin. The parameters of the model are: 
the nucleon mass, of the order of M — 939 MeV, de- 
pending on the parametrization used, the masses of the 
mesons m s , m v , m pi also model dependent, the electro- 
magnetic coupling constant e = \/47r/137 and the cou- 
pling constants between baryons and mesons 

9sB = %sB 9a, 9vB = %vB 9vi 9pB = X pB 9 pi 

where x s b, x v b and x p b are equal to 1 for the nuclc- 
ons and deltas and acquire different values in differ- 
ent parametrizations for the other baryons. In the 
present work we have used x s b = 0.7, x v b = 0.783 
and x p B = 0.783 for couplings between mesons and 
A, £,£, E*,5* and x s b — x vB — x pB — between 
mesons and the fl. 

Non-linear a terms are also included in all parametriza- 
tions through the constants k and A and a non-linear u> 
term is present in the TM1 0] parametrization through 
the constant £. A density dependence is introduced 
through the non-linear a — p and w — p couplings [HI], 
just present in the NLcup model. We have followed the 
prescription of [l7| , where the starting point was the NL3 
parametrization and the g p coupling was adjusted for 
each value of the coupling A$ studied in such a way that 
for hp = 1.15 fm _1 (not the saturation point) the sym- 
metry energy is 25.68 MeV. In the present work we set 
A s = as in [l3j]. Other possibilities for this model with 
(7 — p and u! — p couplings have already been discussed in 
the literature. All coupling constants are adjusted in or- 
der to reproduce the nuclear matter saturation properties 
given in Table HI 



B. Density dependent hadronic model 

The Lagrangian density of the density dependent 
model we use next reads: 

£ = J2$ B W - T vbV» - Y pB r ■ h» 



, e il+jZgg) A p j _ (M 



V'l 



{d^cj) - m 2 ^ 



-mtV u V^--Bp lJ -B 



\™ 2 p h p -W- 



\f,uF^ (2) 



where B M „ and F^ are defined after eq.JTJ). The 

parameters of the model are again the masses and the 
couplings, which are now density dependent, i.e., T s re- 
places g s , T v replaces g v and T p replaces g p . Once again 
the relations 

r.ss = x sB r s , r vB = x vB r.„, r pB = x p b r p 

hold and these density dependent couplings r s , r„ and 
T p are adjusted in order to reproduce some of the nu- 
clear matter bulk properties shown in Table HI using the 
following parametrization: 



Ti(p) = Ti(p sa t)hi(x), x = p/p sat , 



with 



and 



hi(x) 



l + Ciix + di) 2 ' 



S, V 



h p (x) = cxp[-a p (a; - 1)], 



(3) 



(4) 



(5) 



-fl 



with the values of the parameters mi, Ti(p sat ), aj, hi, 
Ci and di, i = s,v,p given in [3|. This model does not 
include self-interaction terms for the meson fields (i.e. 
k = 0, A = and £ = ) as in NL3 or TM1 parametriza- 
tions for the NLWM. 

Two parametrizations of this density dependent 
hadronic models are discussed next. The original one, 
that we refer to as TW and a more recent parametriza- 
tion known as DDME1 also shown to provide a good 
description for the properties of many stable nuclei. 



C. Calculations 

We start by applying the Euler-Lagrange equations 
to the Lagrangian densities and obtaining the mesonic 
equations of motion. Then a mean field approximation 
is enforced and the equations of motion, which have to 
be solved in a self-consistent way become dependent of 
the baryonic densities. At this point the electromagnetic 
field is neglected. For the different parametrizations of 
the non-linear Walecka model, the densities read: 

/A 

with M B = Mb — g S B <t>, B± stands respectively for 
baryons and anti-baryons, E*(p) = \/p 2 + M* 2 and 

f B± = l/{l + exp[(E*(p)Tv B )/T}}, (7) 

where the effective chemical potential is 

vb = Vb - 9vbV - g P B t 3B b . (8) 
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The baryon chemical potentials are computed in terms of 
their quark constituints, i.e., fi n = p, u + 2pj_ and related 
expressions for the others. As particles and anti-particles 
have to be computed separately for the calculation of 
their yields, after the self-consistent calculation they are 
calculated as 



p B+ = (2J B + 1) 



p B - = (2J B + 1) 



d 3 p 



h 



!b- 



(9) 



(10) 



where Jb = 1/2 and 3/2 respectively for the baryonic 
octet and decuplet. For the density dependent hadronic 
model, the expressions are very similar except for the ef- 
fective mass M B = Mb—T s b <j) and the effective chemical 
potential 

vb = Pb — FvbVq — T p b t 3 b bo — T,q TW , (11) 
where the rearrangement term Ti^ TW is given by 



yR 

^0 TW 



E 



dp 



PbVq 



or 



pB 



dp 



i~3b Pb b 



c)T 



sB 



dp 



-PsBPO, 



(12) 



Moreover, as we are interested in obtaining also the 
production of pions and kaons, they are introduced 
through Bosc-Einstcin distribution functions 



2J M + 1 
2tt 2 



p dp 



exp[(Ei - m)/T\ - 1 



(13) 



where i = n + , ir~ , ir°, K + , K~, K°, K°, and the corre- 
sponding vector mesons p and K* , with Jm — and 1. 
Ei = \J p 2 + mf and the chemical potencials are again 
written in terms of their quark constituints, namely, 
fi„+ = Pu — Pd, mu K + — p u — Ps and so on. We have con- 
sidered that they behave like a free gas and their proper- 
ties are not changed due to their interaction with matter 
and, therefore, the fraction of produced mesons is deter- 
mined statistically from their free space properties. 

In order to obtain the particle yields and respective 
densities three conserved quantities are considered: the 
total strangeness is set to zero, the total number of 
baryons in a Au+Au collision is N B = 2{N + Z) = 394 
and the total isospin is I3 = (Z — N)/2 = —39. Our 
code deals with 6 unknows, the three meson fields and 
the three independent quark chemical potentials (p q , q = 
u, d, s), solved in a self-consistent manner. 

Finally, as a last test, we have disregarded the inter- 
actions by setting all couplings equal to zero. We have 
also performed the same test, i.e., for free Fermi parti- 
cles and including the ten baryons with total spin 3/2, 
as suggested in (48[. We next analyse our results. 



III. RESULTS 

As already mentioned, in obtaining our results we have 
fixed the relations between the baryons and the mesons 
as x s b = 0.7, x v b — 0.783 and x p b = 0.783, except 
for the nucleons and deltas, when they are 1 and for the 
f2's, when they are zero. We have verified that other 
choices give slightly different results but have decided to 
fix just one possibility to restrict our test to the model 
differences. 

As stated in the Introduction, we have implemented a 
X 2 fit as in [43j in order to obtain the temperature and 
chemical potential for each model. In Table [TTl we show 
our results for the different models studied corresponding 
to the temperature and chemical potential that for each 
model give the minimum value of the quadratic deviation 
X 2 : 



theo\2 



07 



and 



2 _ v (kt p - n\ 



theo\2 



(14) 



(15) 



where TZ^ P and JZl heo are the i th particle ratio given 
experimentally and calculated with our models and Oi 
represents the errors in the experimental data points. 

In obtaining the best fit values for the temperature 
and chemical potentials, we have used the experimental 
ratios appearing in Table |n] four times for p/p, twice 
for 7r~ /n + and four times for K~ /K + , all with the same 
weight. We have also taken into account the K°* /h~ and 
K°*/h~ ratios, where h~ is the net sum of all negative 
electrically charged hadrons. Instead we could have taken 
the mean value of the measured values and a statistical 
average value of the errors. We have checked that the 
results were similar. 

One can observe from Table [TTJ that we have found in- 
credibly good agreements between experimental data and 
all tested models. The only exceptions are the free octet 
and free octet+decuplet models that cannot describe the 
P/tt~ ratio. Although the inclusion of the decuplet im- 
proves the result, it remains far behind the experimental 
value. These results show that a proper treatment of 
the interaction has to be considered. In [43|], to account 
for the repulsive interaction between particles an eigen- 
volume was assigned to all particles. We have obtained 
approximately the same freeze-out temperature for all 
models, namely, ~ 146.5 — 152.5 MeV. The higher limit is 
close to the value obtained with two parametrizations in 
[48| . As for the chemical potentials, our results are very 
model dependent. While all four parametrizatins of the 
NLWM pb predict a similar chemical potential ~ 47 — 48 
MeV, equal or slighlty lower than the prediction of [48| . 
where the encountered value lies between 48.3 and 54.6 
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MeV, the density dependent models predict larger val- 
ues, respectively 59 and 62.5 for DDME1 and TW. To 
test the importance of including the interaction between 
particles we have set the couplings to zero. We have 
obtained a worse x 2 value compared with the values ob- 
tained when interactions arc included. Although we get 
a similar freeze-out temperature, the chemical potentials 
obtained with no interation are 15-20 MeV smaller. 

In Table II we include other ratios not yet measured. 
We conclude that they are not very sentive to the model 
except fj/fi. The results that are worse reproduced in- 
volve mesons, K — /n— . This seems to indicate that our 
treatment of the mesons, with constituent masses equal 
to the vacuum values, is too naive and an improvement 
should be considered. We have checked that, if we renor- 
malize the mass of the mesons according to the results 
given in [48j namely, consider the pion and kaon mass 1.1 
times the vacuum mass and the K* mass 0.9 of the vac- 
cuum mass we would get \ 2 = 5-3 and T — 149.5 MeV, 
Hb = 48.75MeV for the NL3 parametrization. 

For the sake of completeness, in Table HTT1 we display a 
sample of hadron effective masses. In this table N stands 
for nucleons. One can see that, with very few exceptions, 
they are basically model independent and the influence 
of the medium is rather low for all baryons, excluding the 
nucleons. For nucleons the effective mass is at most 15% 
smaller than the vacuum values. For the other hyperons 
we get 5-6% reduction. These masses seem to show that 
the freeze out occurs below the critical temperature. This 
means that the freeze-out occurs at a temperature below 
the phase transition to a massless baryon phase. 

IV. CONCLUSIONS 

In the present work we discussed the necessity of find- 
ing constraints able to discard or select appropriate rela- 
tivistic models. As all models are parametrized to fix nu- 
clear matter saturation properties, once they are extrap- 
olated to lower or higher densities or finite temperature, 
all sorts of results can turn up. Concerning astronomical 
measurements, observations that could help in ruling out 



possible mass-radius relations are most welcome. Con- 
straints can also come from the measurement of a precise 
neutron skin thickness or from heavy ion collisions, the 
topic of our calculations here. 

By analysing the particle densities and production 
yields obtained in a Au+Au collision at RHIC we have 
verified that all models describe the data in a similar way 
and therefore these data can not be used to establish a 
constraint. This is due to the fact that the chemical po- 
tential and consequently, the baryonic density involved 
in the hadronization process are very low. We have also 
confirmed that in order to improve the fit we must go 
beyond the the naive way mesons were included. 

The analysis of the results of the Pb+Pb collision at 
SPS may be more adequate in the search for constraints 
to different models since, at least, with thermal models, 
the chemical potential seems to be much larger. This 
work is under investigation and, if proved to be an ad- 
equate choice, some improvements can be done for the 
achievement of a better description. The inclusion of 
strange mesons as mediators of the hyperons |53j, re- 
sponsible for explaining the strongly attractive hyperon- 
hyperon interaction as observed in double A hypernuclei 
will certainly influence the densities of particles contain- 
ing strangeness. Also, the inclusion of the scalar isovector 
virtual <5(ao(980)) field, that introduces in the isovector 
channel the structure of relativistic interactions and af- 
fects the behavior of the system at high densities or high 
temperatures [29[ can also bring some modifications. Still 
other relativistic models as the quark-meson-coupling 
model (QMC) [Hi| that presents a quarkionic structure 
inside the hadrons should be investigated. These calcu- 
lations are already under way. On the other hand, more 
data is also required. 
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TABLE I: Nuclear matter properties. 



ratio 


NL3 


TM1 


GM1 


GM3 




NLujp 




TW 


DDME1 




\7] 


m 


m 


m 




ri3] 




m 


® 












A„ = 0.01 A„ = 0.02 


A* = 0.025 






B/A (MeV) 


16.3 


16.3 


16.3 


16.3 


16.3 


16.3 


16.3 


16.3 


16.2 


pa (fm- 3 ) 


0.148 


0.145 


0.153 


0.153 


0.148 


0.148 


0.148 


0.153 


0.152 


K (MeV) 


271 


281 


300 


240 


271 


271 


271 


240 


244.5 


Ssym. (MeV) 


37.4 


36.9 


32.5 


32.5 


34.9 


33.1 


32.3 


32.0 


33.1 


M*/M 


0.60 


0.63 


0.70 


0.78 


0.60 


0.60 


0.60 


0.56 


0.578 
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TABLE II: Comparison of experimental particle ratios and 
relativistic mean field models. 



ratio 


exp. data 


exp 


NL3 


TM1 


GM1 


GM3 


TW 


DDME1 


octet 


octet + 








12 


[8] 


[9] 


[9] 


[3J 


[5] 




decuplet 


P/p 


0.65±0.07 


STAR 


0.650 


0.646 


0.626 


0.597 


0.656 


0.663 


0.661 


0.649 




0.64±0.07 


PHENIX 




















0.60±0.07 


PHOBOS 




















0.64±0.07 


BRAHMS 


















p/n~ 


0.08±0.0f 


STAR 


0.075 


0.072 


0.072 


0.063 


0.076 


0.074 


0.039 


0.041 


1T~ /lY + 


f.00±0.02 


PHOBOS 


0.998 


0.991 


0.999 


1.00 


1.01 


1.01 


1.00 


1.01 




0.95±0.06 


BRAHMS 


















K-/K+ 


0.88±0.05 


STAR 


0.912 


0.911 


0.907 


0.905 


0.896 


0.900 


0.961 


0.941 




0.78±0.f3 


PHENIX 




















0.91 ±0.09 


PHOBOS 




















0.89±0.07 


BRAHMS 


















K /IT 


0.149±0.02 


STAR 


0.234 


0.234 


0.242 


0.243 


0.228 


0.227 


0.232 


0.235 


A/A 


0.77±0.07 


STAR 


0.681 


0.680 


0.666 


0.644 


0.663 


0.675 


0.687 


0.689 


S /S 


0.82±0.08 


STAR 


0.746 


0.747 


0.735 


0.711 


0.739 


0.748 


0.714 


0.732 


K /h 


0.06 ± 0.017 


STAR 


0.058 


0.059 


0.063 


0.064 


0.064 


0.063 


0.060 


0.061 


K°* jhr 


0.058 ± 0.017 


STAR 


0.053 


0.054 


0.057 


0.058 


0.056 


0.056 


0.057 


0.057 


Q/n 






0.693 


0.699 


0.715 


0.723 


0.585 


0.586 


- 


0.784 


Q./-K- 






0.001 


0.001 


0.002 


0.002 


0.001 


0.001 


- 


0.001 


A/h~ 






0.021 


0.020 


0.021 


0.020 


0.023 


0.022 


0.013 


0.014 


Q/E- 






0.172 


0.178 


0.195 


0.211 


0.173 


0.174 


- 


0.250 


\/K°* 






0.351 


0.341 


0.331 


0.301 


0.363 


0.353 


0.226 


0.228 


H-/A 






0.226 


0.226 


0.227 0.221 


0.296 


0.295 


0.241 


0.222 


S-/A 






0.332 


0.332 


0.341 


0.343 


0.330 


0.327 


0.312 


0.322 


S-/K- 






0.047 0.045 


0.050 


0.049 


0.048 


0.045 


0.027 


0.030 


E-/K- 






0.035 


0.034 


0.035 


0.033 






0.019 


0.022 


T(MeV) 






149 


149 


152. 


152.8 


146.6 


146.2 


146.4 


148.8 


Mb (MeV) 






47.5 


46.5 


47.5 


48. 


62.8 


57. 


30.5 


32.5 


p xl0~ 3 (fm~ 3 ) 






8.37 


8.03 


9.62 


9.95 


4.90 


4.45 


2.41 


4.77 


x 2 






23.94 


24.43 


27.99 


33.19 


22.18 


21.83 


45.44 


41.63 


radius (fm) 






22.4 


22.7 


21.38 


21.14 


26.77 


27.65 


33.91 


27.02 



TABLE III: Effective masses obtained with the temperature 
and chemical potentials given in Table ITT1 





NL3 TM1 


GM1 


GM3 TW DDME1 




m 


[8] 


[9] 


\9] 


m 




M* N /M N 


0.88 


0.89 


0.90 


0.93 


0.87 


0.87 


MX/ Ma 


0.93 


0.93 


0.94 


0.96 


0.92 


0.92 




0.94 


0.94 


0.95 


0.96 


0.93 


0.93 


Mi /M B 


0.94 


0.95 


0.95 


0.96 


0.93 


0.94 


Ml /M a 


0.91 


0.92 


0.93 


0.95 


0.90 


0.90 


Af£./Afe. 


0.94 


0.95 


0.95 


0.95 


0.94 


0.94 




0.95 


0.95 


0.96 


0.97 


0.94 


0.94 


M s \/M n 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 



